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ABSTRACT 

We report 4.5 /mi luminosities for 27 nearby (D<5 Mpc) dwarf irregular galaxies measured with the Spitzer 
Infrared Array Camera. We have constructed the 4.5 /im luminosity-metallicity (L-Z) relation for 25 dwarf 
galaxies with secure distance and interstellar medium oxygen abundance measurements. The 4.5 /im L-Z 
relation is 12+log(0/H) = (5.78 ± 0.21) + (-0.122 ± 0.012) M [4 . 5] , where M [4 . 5] is the absolute magnitude at 
4.5 /im. The dispersion in the near-infrared L-Z relation is smaller than the corresponding dispersion in the 
optical L-Z relation. The subsequently derived stellar mass-metallicity (M*-Z) relation is 12+log(0/H) = (5.65 
± 0.23) + (0.298 ± 0.030) log M*, and extends the SDSS M*-Z relation to lower mass by about 2.5 dex. We 
find that the dispersion in the M*-Z relation is similar over five orders of magnitude in stellar mass, and that 
the relationship between stellar mass and interstellar medium metallicity is similarly tight from high-mass to 
low-mass systems. We find a larger scatter at low mass in the relation between effective yield and total baryonic 
mass. In fact, there are a few dwarf galaxies with large yields, which is difficult to explain if galactic winds 
are ubiquitous in dwarf galaxies. The low scatter in the L-Z and M*-Z relationships are difficult to understand 
if galactic superwinds or blowout are responsible for the low metallicities at low mass or luminosity. Naively, 
one would expect an ever increasing scatter at lower masses, which is not observed. 
Subject headings: galaxies: dwarf — galaxies: evolution — galaxies: irregular — infrared: galaxies 



1. INTRODUCTION 

The luminosity-metallicity (L-Z) relation for gas-rich star- 
forming galaxies is a robust relationship which holds over ten 
magn i tudes in galaxy optica l luminosity (e.g., IZaritskv et alJ 
119941 iTremontietalJ 120041) . The physical basis for this 
relation has generally been attributed to a mass-metallicity 
relation for galaxies. The L-Z relation has also become an 
important tool to examine the time-evolution of past chemical 
enrichment and the (stellar) mass-metallicity (M„-Z) relation 
for more massive galaxies at distant epochs, where emis- 
sion lines at rest-frame optical wav elengths are r edshifted 
into the near-infrared; se e, e.g., Kobulnickv & Zaritskv 
(119991): iKobulnickv et all (120031): iLiang et al l (12004 ) ; 
[ Kobulnickv & Kewlevl (12004: IShaplev et alf (12* 004 ) ; 
ICiallazzi et alJ d2005|): iHovosetalJ d2005l): lMaieretal. 1 
l|2005[) : ISavaglio et afld2p05l):IShaplev et al J (120051) :lErb et al.1 
(120061) : lLamareille et alJ ( 120061) : ILiang et alJ (I2006F The 
growing body of observations at intermediate- and high- 
redshifts has led to work on a variety of galaxy formation 
models and their su bsequent predictions for L-Z and M^-Z 
relations (e.g., Somerville & Primack 1999; Bell & Bower 
| 2000i [Tamura e t alJ20011lE>ekel & Wool2003UDe Lucia et al.1 
2004; Savagli o et all2005UTissera et all2005l) . 

For nearby dwarf irregular and other star-forming dwarf 
galaxies, the corresponding L-Z relat ion has tr aditionally 
been studied at optical wavelengths (e.g.J Leaueux et alJl979t 
Skillm anet al.lll989t [Richer & McCallll 19951 ISkillman et alJ 
1997 1 : iPilvuginl 120011 iGarnettl [20021 Iffidalgo-Gamez et aLl 
20031 IH. Leeet alJ l2003albld iPilyugin et alJ 120041 
van Zee & Havnes 2006; vanZeeetal. 2006). If sys- 
tematic uncertainties are minimized, the resulting residuals 
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in the L-Z relation c an be used to examin e effects of 
galaxy evolution (e.g., Kob ulnickv et al.l 120031) . However, 
dispersion in the optical L-Z relation is thought to arise 
from variations in the stellar mass-to-light ratios, which are 
caused by variations in the current star formation rate (SFR) 
amon g galaxies (see, e.g., iBell & de .Tongll2b011 iTumlinsonl 
2006). To minimize the effects of these variations, the L-Z 
relation can be determined at near-infrared (NIR; 1-5 /im) 
wavelengths, where the dominant emission arises from the 
stellar populations. 

Recent ground-based efforts to determine L-Z and M*-Z 
relations for low-mass system s in the NIR have been con - 
str ucted and descr i bed by, e.g.. [Perez-Gonzalez et al.| ([20,0,31), 

e xe et al.1 (120041). ISalzer et aUfcOOSUSayiane et alJ 12005 ), 
veira et alJ (120061) . and Rosenbe rg et al.l ( 120061) . The sen- 
sitivity of the Spitzer Space Telescope provides an excellent 
opportunity to observe total emission from the stellar pop- 
ulations in nearby dwarf galaxies, because most of these 
systems are of intrinsically low surface brightness, and are 
therefore highly challen ging to obs erve in the NIR from the 
ground (e.g., Vaduvescu et al. 2005). The 4.5 /im luminosi- 
ties a re dominated by emission from stellar p hotospheres 
(e.g., iLu et alJ 120031: iDale et alJl2005l lGallianoll2005l) . and 
are less sensitive to bursts of star formation. Observations at 
4.5 /im were chosen also because the passband was designed 
to avoid em ission from polycyclic aromatic hyd rocarbons 
(PAHs; e.g.. iLeger & Pugellll984l lAllamandola et alJll985l) . 
PAH emission does not appear to be sign ificant within lower- 
luminosity metal-poor galaxies (e.g.. lHoucketal.1 120041 
Engelbrac ht et alJl2005l iHo gg et al.ll2005l) . Although emis- 
sion from hot dust might be significant in disk galaxies at 
4.5 /im (e.g., |Roussel et al. 2005), our 8 /im observations 
have shown that the contributions from PAHs and warm dust 
are not significant in our sample of low-mass dwarf galaxies 
( Jacks on et all20 06). 

In light of the M*-Z relation for massive galaxies deter- 
mined by Trem onti et al.l (120041) . we consider nearby dwarf 
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galaxies, which extend the mass range down by roughly 
2.5 dex. This allows us to examine the M»-Z relation over 
five orders of magnitude in stellar mass. To construct the 
L-Z and M*-Z relations at 4.5 /im, we have chosen nearby 
star-forming dwarf irregular galaxies that have distances mea- 
sured from stellar constituents, and oxygen abundances de- 
rived from high-quality spectra (namely, [O III] A4363 mea- 
surements). We have assembled a sample of 27 dwarf galax- 
ies at distances below about 5 Mpc; some properties of these 
dwarf galaxies are listed in Tabled 

The remainder of this paper proceeds as follows. Those 
who wish to skip the details should go directly to § where 
we interpret our key results in the context of the findings by 
iTremonti et al.l J2004I) . Observations and reductions of the 
data are presented in § [2] We describe the measurements 
of flux densities and absolute magnitudes of dwarf galaxies 
at 4.5 /im, and the corresponding NIR L-Z and M»-Z re- 
lations in §§|3]and|4] We have also gathered from the lit- 
erature other L-Z and M*-Z relations for local and distant 
star-forming galaxies. As all of the galaxies in our sample 
have direct measurements of their H I gas, we also examine 
relations with gas-to-stellar mass ratio and H I rotation ve- 
locity in §§|5]and|S] respectively. We give additional dis- 
cussion and a summary in §§0and|8] respectively. For the 
present discussion, we use the notation [O/H] = log(0/H) - 
log(O/H) , where the solar value of the o xygen abundance 
is 12+log(0/H) = 8.66 (lAsplund et al.l2004tlMeIendezl2004l) . 
We adopt Z Q = 0.01 26 as the revised s olar mass fraction in 
the form of metals ( Asnl und et all20 04). 

2. OBSERVATIONS AND REDUCTIONS 

Data were acqu i red w ith the Infrared Array Camera 
(IRAC; Fazio et al. 2004) on the Spitzer Space Telescope 
( Wern er et al.l 120041) during the first year of nominal opera- 
tions from 2003 December to 2004 December. In Guaranteed 
Time Observation (GTO) program 128 (P.I. R.D. Gehrz), ten 
galaxies were imaged in channels 2 (4.5 /im) and 4 (8 /im). 
The Astronomical Observation Request Key (AORKEY) for 
each galaxy is listed in Tabled We performed five dither po- 
sitions per pointing with exposure times of 200 seconds each, 
except for NGC 55 and WLM which had exposure times of 
100 seconds per dither position. A single 5' x 5' pointing was 
used for all targets, except for IC 1613, NGC 55, NGC 3109, 
and WLM, which were mosaicked in array coordinates using 
a 2 x 2 (10' x 10'), a 1 x 2 (5' x 10'), and a 1 x 3 (5' x 15') 
pointing grid, respectively. NGC 55 was mosaicked in sky 
coordinates, resulting in a 15' x 10' field-of-view. Using the 
MOPEX 3 reduction package, the final mosaics were created 
as described in Jack son et alJ f2006). These images contained 
the appropriate astrometry in Epoch J2000 coordinates and 
were flux-calibrated in units of MJy sr" 1 . Postage-stamp im- 
ages at 4.5 /im of ten galaxies from program GTO 128 are 
shown in Fig.[fl note that the galaxies are displayed with the 
same isophotal scaling to emphasize the large range in surface 
brightness. 

Total photometry was obtained with the application of el- 
liptical apertures, and adjusting the dynamic range to show 
the spatial extent of the emission at 4.5 /zm. We also re- 
trieve d optical images fro m NED and the Local Group Sur- 
vey 4 (iMassev et alJl2002l) . and we "stretched" these images 

3 MOPEX is available from the Spitzer Science Center at 
http://ssc.spitzer.caltech.edu/postbcd/ 

* More about the Local Group Survey is found at 



visually (in intensity scale) in order to discern the outer limit 
of the optical extent and to compare the spatial extents of the 
optical and NIR luminosities. Elliptical apertures were con- 
structed to optimize the amount of light measured for each 
target galaxy. Six to seven circular apertures in the immedi- 
ate vicinity of the elliptical aperture were used to sample and 
determine the mean background level surrounding the galaxy. 
This method benefits from the sampling of both the sky and 
the foreground stellar contamination. An example is shown in 
Fig. |2 We used the IRAF 5 task polyphot to perform the aper- 
ture photometry and measure flux densities. A variety of tests 
were performed to mask and remove negative pixels, bright 
foreground stars, and bright pixels caused by "bleeding" from 
bright stars. Additional experiments (varying the size of the 
elliptical aperture; increasing the number of circular "back- 
ground" apertures) were used to test our final photometry re- 
sults. The subsequent range in derived luminosities falls well 
within the quoted uncertainties for the distance modulus and 
the total errors for IRAC flux measurements (see below). 

We have also added from the Spitzer public data archive 
eleven dwarf galaxies from the SINGS Legacy program 
( Kennicutt et alJ l2003bl) . one dwarf galaxy from program 
GTO 59 (P.I. G. Rieke), and five dwarf galaxies from program 
GTO 69 (P.I. G. Fazio). To ensure homogeneity in the derived 
absolute magnitudes at 4.5 /im, we have performed the same 
procedure described above to determine total fluxes and lu- 
minosities. As shown in Fig. [5] our measured flux densities 
for the se lected SINGS ga laxies are consistent with those re- 
ported bv lDale et ail J2005I) . 

We adopt 10% as the error in the total flux, assuming this 
is the maxim um error of the absolute calibration in channel 2 
dHora et al.l2004l) . Various image-based corrections (array lo- 
cation dependent photometric correction, pixel-phase correc- 
tion) in channel 2 are well below th e uncertainty in the ab- 
solute calibration jReach et al.l20 05). and are not included in 
this paper. We have applied a correction to account for ex- 
tended surface brightness by multiplying measured channel 2 
fluxes by a factor of 0.937 dReach et aDl2005l) to obtain the 
appropriate scaled results. With our measured flux densities 
at 4.5 /im, we derived absolute magnitudes, M[4.5], by using 
the adopted zero-mag nitude flux density, = 179.7 ±2.6 
Jy (Reach et al. 2005). We have also assumed zero extinction 
at 4.5 urn (e.g.. lLutzH99'9tlIndebetouw et alJ2005l) . and neg- 
ligible contribution by diffuse H I Brack ett-a em ission to the 
overall galaxy flux (e.g.. lLu et al.l2003l 020041) . 

A comparison between B and 4.5 /im luminosities is shown 
in Fig. |4] where B luminositie s corrected for G alactic extinc- 
tion are taken from Kar achentsev et all J2004I) . For the re- 
mainder of this paper, we will refer to IRAC channel 2 lu- 
minosities as "[4.5]". Galaxy luminosities at 4.5 /im are also 
plotted against B-[4.5] color; no correlation is evident. The 
mean B-[4.5] color for our sample in Tablets 2.08 ±0.59, 
which is consistent with the range of observed B-K colors 
for a sample of ne arby late-type dwarf galaxies in 2MASS 
dJarrettetal 120031). 

3. THE NEAR-INFRARED LUMINOSITY-METALLICITY 
RELATION 

|http : / / www . lowell ■ edu/users/massey/lqsurvey ■ html| 

5 IRAF is distributed by the National Optical Astronomical Observatories, 
which are operated by the Associated Universities for Research in Astron- 
omy, Inc., under cooperative agreement with the National Science Founda- 
tion. 
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Twenty-five dwarf galaxies in our sample have reported 
oxygen abundances, and 21 dwarf galaxies have oxygen 
abundances derived from [O III] A4363 measurements (Ta- 
ble 0. In the absence of [O III] A4363, oxygen abun- 
dances for DD O 53, Ho I, an d NGC 3109 were de- 
rived using the McGaugh| d!991l) bright-line calibration 6 . 
We adopt for Peg DI G the oxygen abundance derived by 
ISkillman et aD (l99% . We have assumed oxygen abun- 
dances are spatially homogene ous in dwarf irregular g alaxies 
(e g , [K obulnickv & Skill manll 19971 iKobulnickv et alJll997t 
IH. Lee & Skillmarl2004tlH. Lee et all2005U2006al) . 

To obtain fits between absolute magnitude and oxygen 
abundances for our sample of dwarf galaxies in B and 
at 4.5 /mi, we used the ordinary least squares bisec- 
tor method (see Akritas & Bershadv 1996, and references 
therein). Specifically, these are bivariate least-squares fits 
with the bisector method incorporating heteroscedastic errors 
(in both luminosity and oxygen abundance). The two fits are 
expressed as 

1 2 + log(0/H) = (5 .94 ± 0.27) + (-0. 1 28 ± 0.0 1 7) M B ( 1 ) 
= (5.78 ± 0.21) + (-0.122 ± 0.012)M [45] ,(2) 

where Mb is t he absol ute magnitude in B, corrected for Galac- 
tic extinction ( Karachentsev et al. 2004), and M\\ 5] is the ab- 
solute magnitude at 4.5 /im, respectively. The fits in B and 
[4.5] are shown in Fig. [5] The slopes, intercepts, and their 
corresponding errors are computed with the bi-variate corre- 
lated errors and scatter regression code 7 , where we have as- 
sumed uncorrected errors and set the covariance to zero. The 
resulting errors incorporate uncertainties associated with the 
data points being fitted. All of the fits below use the same 
technique, and we refer to this procedure as the "bisector" 
method 8 . The fit at B is con s istent with similar f i ts ob tained 
bv.e.g.JH. Lee et al] J2003bl) . lvan Zee & Havnesl J2006I) . and 
van Zee et alJ d2006). The NIR slope is marginally smaller 
than the optical slope, although the derived slopes and inter- 
cepts for both L-Z relations are consistent with each other. 

Variations in the stellar mass-to-light ratio can be responsi- 
ble for some or all of the scatter in the L-Z relation. Moreover, 
the optical luminosity of a fainter galaxy is affected more by 
a burst of new star formation t han the luminosity of a brighter 
galaxy. Skillm an et al.l {1997) proposed a comparison of the 
residuals in luminosity against color to test whether variations 
in the stellar mass-to-light ratio were the main cause of the 
scatter. A limited set of data at the time were plotted against 
optical B-V color, and were suggestive of a correlation be- 
tween the residuals in B luminosity versus B-V. 

We show now that the scatter in the L-Z relation dimin- 
ishes when luminosities are measured at longer wavelengths. 
Residuals in oxygen abundance are plotted against B and [4.5] 
absolute magnitude in the top and bottom panels, respectively, 
of Fig. [5] The dispersions in the optical and NIR L-Z relations 
are 0.161 dex and 0.122 dex, respectively. The scatter in the 
NIR L-Z relation is reduced compared to the optical L-Z re- 
lation. The dispersion in the L-Z relation was expected to 
decrease at longer wavelengths, as NIR luminosities are less 

6 A nalytical expressions for the calibration are found in Kobulnickv et al. 
1 1999). Measurements of [N II]/[0 II] (< 0.1) break the double-valued de- 
generacy and point to the "lower bran ch," which is appropriate for dwarf 
galaxies; see H. Lee & Skillman 1 2004) for additional discussion. 

7 The code is available from STATCODES: 
|http : / / www . astrostatistics .psu . edu/stat codes /| 

8 For each fit, we perform 10000 "bootstrap" simulations to converge onto 
a solution with the appropriate errors. 



sensitive to extinction by dust and variations in the present- 
day SFR. 

The smaller slope and the reduced dispersion in the NIR are 
in agreement wit h the results for sta r-forming galaxies in the 
field reported by Sa lzer et all (I2005I) . Higher-luminosity and 
highe r-mass gala xies contain more metals and more dust (e.g., 
iRosenberg et al.l2006l) : these galaxies would be optically un- 
derluminous, which would yield a steeper slope in the optical 
L-Z. The effects of absorption by dust vary from galaxy to 
galaxy, which introduces greater overall scatter in the optical 
L-Z. Absorption by dust is mitigated in the NIR, which pro- 
duces an L-Z relation with a shallower slope and less scatter. 
For our present sample of dwarf galaxies, we have assumed 
that dust is less important (see §0 than the variations in stel- 
lar mass-to-light ratios. 

We now examine residuals in oxygen abundance and lu- 
minosity against B-[4.5] color. The results are shown in 
Fig El Residuals in oxygen abundance from the B L-Z rela- 
tion are correlated with color in the sense that redder galaxies 
have larger abundances compared to the best fit. On the other 
hand, there is very little to no trend with color for the residuals 
from the [4.5] L-Z relation. Residuals in magnitude from the 
optical and infrared L-Z relations exhibit similar behavior as 
above; redder galaxies are fainter in B and no trend is found 
at 4.5 /im. The plots in Fig.[7]clearly show that the dispersion 
is clearly reduced at NIR wavelengths. 

We compare in Table|2]the 4.5 /im L-Z relation with recent 
optical and NIR L-Z relations from the literature. The slopes 
from our optical and NIR relations are smaller than values de- 
rived from similar relations for other samples of nearby dwarf 
galaxies. In addition to listing the slope and intercept val- 
ues for each L-Z relation, we have also computed the oxygen 
abundance from the specified relation at a fiducial galaxy lu- 
minosity equal to M\ = -16. The key result from the compar- 
ison is that our NIR L-Z relation exhibits the smallest disper- 
sion. 

4. STELLAR MASS-METALLICITY RELATION 

The NIR L-Z relation reflects a fundamental relation be- 
tween underlying stellar mass and metallicity, and thus, NIR 
luminosities of dwarf galaxies can d irectly be related to th e 
mass of their stellar populations (e.g., Vaduvescu et al. 2005). 
First, to derive the 4.5 /im luminosity in solar units, we have 
assumed that the absolute magnitude of the Sun is +3.33 in K 
JBell & de J ong 2001), and that the ground-based K-M color 9 
is about for a G2 dwarf star (e.g.. iBessell & Brettl l 1988). 
Thus, we assume that the absolute magnitude of the Sun at 
4.5 /im is ~ +3.3. Even if the NIR luminosity is dominated 
by giant stars, the most extreme Zf-[4.5] colors are between 
+0.1 and -0.1. As we show below, a non-zero /f-[4.5] color 
introduces a zero-point offset in the M*-Z relation. 

We consider an appropriate stellar mass-to-NIR-light ratio. 
We first assumed a constant stellar mass-to-[4.5]-light ratio of 
0.5, as motivated by m odels for the formation of disk galaxies 
(e.g.. lvan den Boschll2002l) . This value of the stellar mass-to- 
light ratio is assumed to be independent of color. The result- 
ing fit to the stellar mass-metallicity relation is 12+log(0/H) 
= (5.26 ±0.27) + (0.332 ±0.033) log M*, and the resulting 
dispersion is 0.125 dex. However, stellar mass-to-light ratios 
can vary by a factor of ~ 2, even at NIR wavelengths (e.g., 
iBell & de.Tong lbOOl). Fortunately, stellar mass-to-NIR light 

9 We assume that the effective wavelength of M is 4.8 /im 1 Cair roins et alJ 
H98l . 
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ratios are not strongly dependent upon the age of the dominant 
stellar population. 

For the next step, we derived stellar mass-to-light ratios 
by accounting for col or variations. We used the models by 
iBell & de Jong ! (1200 It their Tables 3 & 4), where the authors 
developed methods to derive better estimates of stellar mass 
in disk galaxies and to det ermine the baryonic Tu lly -Fisher re- 
lation. For each model in lBell & de Jond (120011) . we derived 
a mass-to-light ratio as a linear function of B-K color, and 
computed the stellar mass for each galaxy. Using the color- 
based functions for B—V and V—K, we applied a simple linear 
combination to obtain the stellar mass-to-light ratio as a linear 
function of B-K color. We write 

log(M*/L K ) = a'+b'(B-K), (3) 

where a' and b 1 are given in Table|3] and B-K = (B— [4.5])- 
(A"-[4.5]) using our measure d fl-T4.51 colors. For a small 
sample of late-type galaxies, Pahr e et all (|2004) determined 
A"-[4.5] in the range between -0.1 and +0.4. We adopt here 
AT-[4.5] = +0.2, and subsequent stellar masses are 0.15 to 
0.16 dex smaller than stellar masses derived with /T— [4.5] = 0. 
We obtain stellar masses with 

log M* = log (M« /Lk) + log L K 

= log(M*/L K )+ [logL [4 .5]-0.4(^-[4.5])] , (4) 

with the stellar mass-to-light ratio defined above. With this 
empirical prescription for stellar masses, we subsequently 
performed a M*-Z fit to each model. The results are shown 
in Table [3] Although stellar masses for a given galaxy can 
vary by as much as ~ 0.5 dex from model-to-model, the sub- 
sequent M»-Z fits spanning a range of about four decades in 
stellar mass are similar for all models. The M*-Z slopes for 
models at subsolar metallicity are marginally smaller than the 
scaled Salpeter stellar initial mass function (IMF) models at 
solar metallicity. 

For this work, we used model 8 in Table [3] (Salpeter IMF, 
subsolar metallicity) and Equations Q and (0} to derive ini- 
tial stellar masses. We then subtracted 0.25 dex from all of 
the stellar masses to scale the masse s from the [S^lnejgr|(|1 955l) 
IM F to thelChabrierl ll2003l) IMF. iTremonti et al.1 (120041) used 
the Kroupa) (120011) IMF to derive stellar masses, which are 
about 0.30 dex smaller than masses derived with the Salpeter 
IMF: lBell et al.l d2003h provide additional discussion on using 
different IMFs to derive stellar masses in late-type galaxies. 
We adopt the errors in [4.5] absolute magnitude as the errors 
in the logarithm of the stellar mass. We note that our derived 
stellar mass for NGC 6822, log (M*/M Q ) = 8.10, is consistent 
with the value (8.24) derived from an independen t model de- 
scribi ng the photometric and chemical evolution (Cari gi et al.l 
2006]). 

The final M*-Z relation is expressed as 

12 + log(0/H) = (5 .65 ± 0.23) + (0.298 ± 0.030) logM* . (5) 

The dispersion is 0.117 dex. In Fig. [8] we show the data and 
the fit expressed by Equation (0 in panel (a). We find an 
excellent correlation between oxygen abundance and stellar 
mass; the Pearson correlation coefficient is r = +0.90. No 
trend is seen in a plot of the residuals in oxygen abundance 
from the M»-Z relation against stellar mass in panel (b). 

We compare our M*-Z relat ion with the relation recently 
derived for massive galaxies. Tremonti et al] d2004l) derived 
the M»-Z relation for over 53000 massive star-forming galax- 
ies at redshifts z < 0.1 from the Sloan Digital Sky Survey 



(SDSS). In Fig. [9] we show the full range in both axes, span- 
ning over two orders of magnitude in oxygen abundance, 
and five orders of magnitude in stellar mass. The data pre- 
sented here extend the SDSS M»-Z relation to lower masses 
by roughly 2.5 dex. The dispersion for our M*-Z relation is 
comparable to that of the SDSS M*-Z relation (w 0. 1 dex). If 
dwarf galaxies have insufficient total masses to retain ejecta 
from supernovae, we expect that the dispersions in both L-Z 
and M*-Z relations should increase with decreasing luminos- 
ity and mass. We will discuss the implications of our L-Z and 
M*-Z results in §□ 

The polynomial fit to the SDSS galaxies is effective in the 
mass range between 10 s 5 and 10 Mq, whereas our rela- 
tion overlaps the SDSS relation in the mass range between 
10 8 - 5 and 10 9 - 5 M Q . At these stellar masses, SDSS galaxies 
are about 0.2-0.3 dex more metal-rich than the present sam- 
ple of dwarf galaxies. We describ e here t he m ost likely rea- 
sons for the offset; see also Savagli o et alJ(l2005l) . First, SDSS 
fibe r spectra pre ferentially sample the inner w 25% of galax- 
ies ( Tremo nti et al.l l2004). Giv en that nearby spiral s exhibit 
radial metallicity gradients (e.g., Zaritsk v et al.ll994l) . derived 
abundances from SDSS spectra might be too high for a num- 
ber of the galaxies. Second, oxygen abundances for the SDSS 
data were derived using bright-line methods in the absence of 
the temperature-sensitive [O III] A4363 emission line. Recent 
work with nearby spiral galaxies has shown that bright-line 
oxygen abundances overestimate the true oxygen abundances 
by roughly 0.2 to 0.3 dex at metallicities solar and above (e.g ., 
Stasihska 2002; Kennicu tTetal]l2003at iBresolin et all 2004). 
These effects together can account for the offset between the 
present sample of dwarf galaxies and the SDSS sample of 
higher-mass galaxies. 

Fig. |2l suggests two intriguing possibilities: (a) over a large 
range in mass (~ 5.5 dex) there is a single and continuous 
M»-Z relation which turns over or flattens at high galaxy 
mass, or (b) the M»-Z relations describing dwarf and giant 
galaxies may have differ ent slopes; this possibility has also 
been s uggested by e.g.,[Salzer eTatl 12005): iRosenberg et alJ 
(2006). It will be very interesting to investigate whether 
slopes of the L-Z and M*-Z relations exhibit a "break" at 
some characteristic luminosity and mass, respectively 10 . 

In Table0] we list the M»-Z relation from the present work, 
as well as other relations from the literature. With increas- 
ing numbers of intermediate- and high-redshift galaxies with 
measured luminosities and abundances, the interpretation of 
the M*-Z relation for massive systems at distant epochs will 
depend upon the specific bright-line calibration used to de- 
rive oxygen abundances. Although selection criteria and the 
methods used to derive masses and metallicities may vary, 
comparisons which begin to tie relations between high-mass 
and low-mass galaxies and between local and distant galaxies 
may help illustrate processes acting on galaxies over a range 
in mass and redshift. 

5. RELATIONS WITH GAS-TO-STELLAR MASS RATIO 

We consider the gas-to-stellar mass ratio, whichjs repre- 
sentative of the baryonic gas fraction (e.g., iGarnettl 120021 

10 Examination of the latest SDSS data (~ 80000 objects in DR4; C. 
Tremonti, personal communication) shows that the subsequent M,-Z fit is 
effective to another 0.5 dex lower in mass to log (M^/Mq) ~ 8. Despite 
the different methods by which stellar masses and oxygen abundances have 
been determined, there appears to be good agreement at log (M*/Mq) ~ 8 
between the upper-end of our M»-Z relation and the lower-end of the SDSS 
M*-Z relation. 
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| H. Lee et al.ll2003bt Islilirnan et alJl2003l Ivan Zee & Havnesl 
2006). These parameters reveal important information about 
the conversion of gas into stars and metals, and about possible 
gas flows. Assuming zero gas flows and instantaneous recy- 
cling, the closed-box model of chemical evolution predicts the 
metallicity by mass Z (Schmidt 1963; Searle & Sargent 1972) 

Z = yln/T 1 , (6) 

where y is the yield and fi = M ?as /(M eas + M») is the fraction 
of total baryonic mass in the form of gas. From the equation 
above, we can derive the oxygen abundance by number 

12 + log(0/H)= 12 + log(y /C) + log(ln [l + (M gas /M»r I ] ) , 

(7) 

where yo is the oxygen yield by mass, £= 11 .728 is the factor 
to convert abundance by mass to abundance by number, M gas 
= 1.36 M m is the total gas mass (which includes helium), and 
M gas /M» is the gas-to-stellar mass ratio, which is independent 
of distance. 

All of the galaxies in the present sample have direct mea- 
sures of their H I gas conte nt. We have adopted 21 -cm flux in- 
tegrals from the catalog by Karachentsev et al. (2004). How- 
ever, we have not scaled the H I mass to the spatial extent of 
the optical or NIR luminosity (re. stellar mass) 11 . The 21-cm 
flux integrals and the derived gas-to-stellar mass fractions are 
listed in Table|5] 

The results are shown in Fig. EH In panel (a), we have plot- 
ted oxygen abundances against the gas-to-stellar mass ratio, 
and the subsequent least-squares fit with the bisector method. 
The fit is shown as a solid line, and is expressed as 

12+log(O/H) = (8.21±0.07)+(-0.435±0.085)log(M gas /M*). 

(8) 

The dispersion is 0.24 dex. In the center panel (b), we have 
plotted outflow models as variations of the closed-box model; 
we assume negligible inflow of metal-poor gas. The pur- 
ple and blue set of c urves, respectively, describe yn = 01 
( Trem onti et al 12004 an d y = 2 x 10"* iU. Lee et alJ2 003b; 
van Zee & Havnes 2006). For each set of curves, the outflow 
rates from right to left are 0, 5, 10, 25, and 50 times the star- 
formation rate. These results suggest reduced yields and/or 
significant outfl ow rates, w hich previous authors (e.g., Garnett 
2002; van Zee & Havnes 2006) have also indicated. 

In Fig. HOt. we have plotted the B-[4.5] color against the 
gas-to-stellar mass ratio, where there is a very clear correla- 
tion. We first point out three outliers in the data. At B-[4.5] 
w +1 and log (M gas /M„) w -0.4, NGC 1569 is unusually 
blue for its gas-to-stellar mass ratio, as this dwarf galaxy is 
known to ha ve just undergone a very st rong burst of star for- 
mation (e.g.. Hsraeill988l tWaiIerlll991l) . At the other end of 
the plot (B-[4.5] w +2, log (M gas /M*) « +1.5) are DDO 154 
and Ho IX. DD O 154 is known to have a very la rge H I halo 
dCarignan & Beaulieulfl989l ICarignan & Purtonlll998l) . Be- 
cause Ho IX is situated (in projection) very close to the spiral 
galaxy M 81, the H I for Ho IX is most likely a part of the 
overall gas distribution for M 8 1 . 

The bisector fit to the data (excluding DDO 165, Ho IX, and 

1 1 The proper interpretation of effective yields with respect to the "inner" 
and "outer" gas relative to the optical extent depends very strongly on the 
assumed metallicity of the outer gas, which is not known or well constrained 
in many nearby dwarf galaxies; however, see, e.g., Skillman et al. 1 1989); 
Kennicutt & Skillman 1 200 1 ); Garnett 1 2002 ) for additional discussion. A 
detailed comparison between "non-scaled" and "scaled" H I content relative 
to the optical extent is outside the present focus of this paper. 



NGC 1569) is shown as a solid line, and is expressed as 

(fi - [4.5])= (2.79 ±0.11) + (-0.89 ±0.13) log(M gas /M»), 

(9) 

and the dispersion in color is 0.32 mag. The fit also implies 
that Lb/L[4.5] oc (M gas /M») ' 36 . The strong correlation may 
simply reflect a more dominant process of the conversion of 
gas into stars. Galaxies with lower gas-to-stellar mass ra- 
tios have already depleted their fuel, and appear redder when 
present-day star formation no longer becomes active. 

We consider possible departures from the closed-box model 
by inverting Equation l|6) to obtain the effective yield, y e g = 
Z/ln(/i _1 ). The effective yield is determined simply from 
the measured oxygen abundances and baryonic gas fractions, 
and is equal to the true yield when the closed-box assump- 
tion applies. Effective yields which are smaller than the true 
yield sugg est ei t her outflow or infall of metal-poor gas (e.g ., 
Edmunds' 1990; iDalcantonl 120061 Ivan Zee & Havnesl 2006). 
We have listed derived values of the effective yield in Table[5] 

In Fig. II 1 h . we plot the effective yield a s a function of 
the total baryonic ma ss, following Fig. 8 in iTremonti et alJ 
(2004): lGaniet3(l2()02l) also obtained a similar plot for a sam- 
ple of nearb y spiral and irregular galaxies. The empirical fit 
obtained by (Tremon ti et al.l l2004) is based upon the premise 
that galactic winds wer e responsible for metals loss; see 
also IVeilleux et alJ l2005 for additional discussion. From the 
present sample, there are dwarf galaxies with effective yields 
about a factor of ten lower than the "asymptotic" yield (log 
Jeff ~ —1.98) for more massive galaxies in the SDSS. We have 
also marked for compari son the solar oxygen mass fr action 
(log Z ,o = -2.27, Aspl und et aDl2004 lMelende z1l2004 and 
the e ffective yield fo r the nearby spiral M51 (log y e g = -2.49, 
iBresolin et alJ20 04) 12 . Accounting for the ~ 0.2-0.3 dex off- 
set described above in §|4] lowers the "asymptotic" yield in 
better agreement with the solar mass fraction and the effective 
yield for M5 1 . We note that these values are also comparable 
to the range of effective yields for dwarf galaxi es (e.g., log 
Veff ~ -2.7 to -2.1) previously r eported by, e.g.. lH. Lee et alJ 
2003bt lvan Zee & Havnesl2006l) . 

Although there may be a trend between yield and baryonic 
mass for the low-mass galaxies, there is however large scat- 
ter, which is about 1.5 dex at ~ 10 8 5 M Q . Three galaxies 
with the largest effective yields (log y & s > -2) in our sam- 
ple are DDO 154, NGC 3109, and Sextans A. However, their 
measured luminosities and abundances are in good agreement 
with other galaxies at comparable abundances and luminosi- 
ti es, respectively. 

Ivan Zee & Havnesl (120061) observed a trend between effec- 
tive yield and H I-gas-to-fi-light ratio for a sample of nearby 
UGC dwarf galaxies, where the H I-gas-to-B-light ratio was a 
measure of the gas-to-stellar mass ratio. These authors ruled 
out infall of metal-po or gas as a mech anism to decrease ef- 
fective yields (see also Dalcanton 2006). In Fig. II lb. we plot 
the effective yield as a function of the gas-to-stellar mass ra- 
tio. Gas-rich systems eventually reach the "asymptotic" yield, 
while gas-poor systems have lower effective yields. The ef- 
fective yield is correlated with the gas-to-stellar mass ratio, 
and the fit is expressed as 

logoff = (-2.75 ± 0.06) + (0.700 ± 0.057) log(M gas /M*) 

(10) 

1 2 For a sample of H II regions in M5 1 , Bresolin et al. 1 2004 ) measured 
electron temperatures from [N II] A5755 measurements, and derived chem- 
ical abundances. The subsequent effective yield for M51 was found to be 
lower by a factor of four than previous estimates. 
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with dispersion in log yield equal to 0.22 dex. In the limit of 
high gas fraction, the closed-box model predicts that the slope 
in equation above is equal to one. 

6. RELATIONS WITH ROTATION VELOCITY 

For a sample of dwarf irregular and spiral galaxies. iGarnettl 
(2002) found a relationship between metallicity and rota- 
tion vel ocity with an apparen t "break" in the relation at 120 
km s" 1 . Pilvugin et al. (2004) peformed a homogeneous anal- 
ysis using his bright-line calibration, and found a similar 
break at 80 km s" 1 . We plot various quantities against H I 
rotational velocity in Fig. 1121 where the range of velocities 
for the present sample is between 10 and 80 km s" 1 ; see Ta- 
ble [2 As an alternative to velocity on a linear scale, we plot 
the logarithm of the velocity. 

In Fig.ll2h. there appears to be a good correlation between 
stellar mass and rotation velocity (r = +0.66), although there 
is a paucity of data at low rotation velocity. The bisector fit is 

logM* = (4.14 ± 0.57) + (2.33 ± 0.36)logv rot , (11) 

and the dispersion in log stellar mass is 0.63 dex. The fit 
can be inverted to give v 10t oc M*° 43 ; the nume rical exponent 
here i s larger than the predicted value of 0.2 ( Dekel & Woo 
2003) by 3.4cr. The observational result suggests either larger 
rotation velocities at a given stellar mass, or smaller stellar 
masses at a given velocity. 

In panels (b) and (c), respectively, we have plotted the gas- 
to-stellar mass ratio and the effective yield against rotation 
velocity. Neither panel reveals a strong trend, although the 
effective yield appears to rise slightly with rotation veloc- 
ity. The rise of th e effective yiel d with rotation veloc ity has 
been discussed b y Garnetn (120021) . IPilvugin et al. (2004|), and 
iDalcantonl (120061) . Galaxies with v ro t ^100 km s -1 have ef- 
fective yields approaching that of the asymptotic yield. Less 
massive galaxies with v ro t < 100 km s" 1 have smaller effective 
yields which are roughly consistent with a solar-type yield, as 
seen above in Fig. The present work suggests that the 
effective yield may be roughly constant, or at the very least, 
is weakly dependent upon rotation velocity. While additional 
data are required for confirmation, this result suggests that 
within the mass range of dwarf galaxies, the effective yield de- 
pends more strongly upon the relative gas fraction than upon 
rotation velocity. 

7. DISCUSSION 

We have shown that the scatter in the optical L-Z relation 
is constant (~ 0.16 dex) over a range of 11 mag in B lumi- 
nosity, and the scatter decreases to 0.12 dex over 8 mag at 
NIR wavelengths. Our subsequent stellar mass-metallicity re- 
lation extends the dynamic range of the SDSS M*-Z relation 
to lower masses by roughly 2.5 dex. We have also shown that 
the scatter for the M»-Z relation appears to be similar over 
5.5 orders of magnitude in stellar mass. 

Models of galaxy formation have been successful in pre- 
dicting the slopes of the L-Z and M*-Z relations. These 
models generally incorporate feedback, which describes how 
effectively massive stars can inject energy int o the surround- 
ing interstellar medium. Dekel & Silk ( 1986) predicted the 
scaling relation, L oc Z 2 7 , for dwarf galaxies formed and em- 
bedded within dark matter haloes, although their model is 
less applicable to dwarf galaxies which have undergone a low 
and constant level of star formation over a Hubble time. The 
derived slopes of the L-Z relations (our Equation |2j are in 



approximate agreement with the lDekel & Silkl (H9 86) predic- 
tion, although our [4.5] L-Z sl ope is smaller t han the predicted 
value by about 2a. Recently, Dek el & Wool ( 120031) predicted 
the scaling relation Z oc M* . The slope of the M*-Z re- 
lation in Equation (0 is smaller th an the predic t ed va lue by 
3.4cr. We note, however, that the lDekel & Wool (120031) scal- 
ing relation was derived for stellar iron metallicities under the 
assumption of constant [Fe/O], which may not be applicable 
in dwarf irregular galaxies (e .g.. Smith et al. 2002). 

Som erville & Primackl (119991) compared predicted L-Z re- 
lations among a number of semi-analytical models. Although 
the details varied with the specific feedback recipe within the 
models, the scatter in the L-Z relations increased with de- 
creasing total luminosity (their Fig. 18). Unfortunately, none 
of the models provided acceptable fits to all of the known ob- 
servational parameters, and these authors stated that the fail- 
ure of the models is likely caused by a strong difference in 
efficiency between ejec ting gas and ejecting metals. 

Models by De Lucia et al. ( 2004) produced M»-Z relations 
consistent with the SDSS M„-Z relation over a simil ar range 
in stell ar mass, down to log (M*/M©) ~ 8.5. iRobertson et all 
(2005) developed models to describe chemical abundance 
patterns observed in the Galactic halo, Local Group dwarf 
spheroidals, and Local Group dwarf irregular galaxies. They 
incorporated chemical enrichment histories and various gas 
accretion histories within cosmologically motivated ACDM 
models for the formation of dwarf galaxies and the Galac- 
tic stellar halo. Th eir models were able to reproduce the 
Dekel & Woo (2003) scaling relation and the observed abun- 
d ance patterns in nearby dwarf galaxies. 

Tremonti et al. (2004) explained that the relationship be- 
tween stellar mass and metallicity was a consequence of the 
relationship between effective yield and total baryonic mass. 
That is, higher-mass systems are able to convert gas into stars 
which produce galaxies with large stellar masses and metallic- 
ities; whereas, lower-mass systems convert less gas into stars, 
and lose their metals as outflows by virtue of the lower to- 
tal masses. These authors state: ". . . It is our view that the 
strong positive correlation between effective yield and bary- 
onic mass is most naturally explained by the increasing po- 
tential barrier that the metal-laden wind must overcome to 
achieve 'blowout.' ..." In fact, since a number of the low- 
mass galaxies in our sample do not depart significantly from 
the L-Z and M»-Z relationships and have effective yields 
close to or exceeding the asymptotic value, it is difficult to ex- 
plain these relationships as a result of very disruptive blowout 
events. 

Present-day galactic winds removing material from low- 
mass galaxies should have dramatic effects on the prop- 
erties of such systems. We consider a prototype dwarf 
starburst gal a xy, NG C 1569, which is undergoing blowout. 
Mart in et al.1 {2002) showed in their Chandra observations 
that most of the oxygen produced in the star burst is in the en- 
riched outflow. Ho wever, the resul t s from llMartin et al.l2002l 
their Table 6) and Anger etti et al.l ( 120051) reveal that at most 
2% of the stellar mass is produced in the burst. Although 
the stellar mass (« 10 8 5 M Q ) of NGC 1569 is at the upper 
end of the range for our sample of dwarf galaxies, an equiv- 
alent NGC 1569-like burst would produce a more significant 
amount of stars in lower-mass systems. However, this high 
level of star formation is not presently seen in the other low- 
mass galaxies from our nearby sample, and cer tainly not see n 
in dwarf irregulars within the Local Group (cf. Mateo 1998). 
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We hypothesize that a less energetic form of metal- 
enhanced mass loss than blowouts could explain the small 
scatter in the L-Z and M*-Z relationships. Perhaps the slow 
loss of metals from an x-ray corona which is fed by individual 
supernova is sufficient to explain the observed relationships. 
Since this mechanism is an average over many events instead 
of depending on a small number of very energetic events, a 
smaller scatter would arise naturally. A low efficiency of star 
formation in dwarf irregular galaxies could account for a gen- 
tler form of mass loss. 

While a number of the more massive dwarf galaxies in the 
present sample h ave experienced or w i ll experience an out- 
flow event (e.g., Summers i ari1l2003 lH artwell et al] l2004l 
iGrimes et alJl2005l lOtt et alJl2005l) . the challenge remains to 
explain the uniform scatter in both the NIR L-Z relation over 
a range of about 1 1 magnitudes in luminosity, and the M*-Z 
relation over a range of about 5.5 decades in stellar mass. We 
suggest the following lines of research for the near future. 
New generations of models for galaxy formation which in- 
corporate varying degrees of galaxy outflows should be able 
to predict the slope and the scatter in the L-Z and M*-Z re- 
lations. Large uniform surveys of galaxies are required to 
examine the percentage of sta rs formed in starbursts as a 
function of galaxy mass (e.g., Perez-Gonzalez et al. 2003). 
Resul ts derived from the resolved stars in NGC 1569 (e.g., 
Anger etti et al.fe 005) are promising, because similar studies 
will allow us to determine how much stellar mass has formed 
in the most recent bursts to a lookback time of a few Gyr. Ad- 
ditional NIR imaging and spectroscopy for a larg e number of 
dwar f galaxies within the Local Volume ( Karachen tsev et all 
120041) can test whether the scatter in the L-Z and M*-Z re- 
lations is truly comparable to the scatter in the respective re- 
lations for massive galaxies. With growing interest in under- 
standing further the stages of galaxy assembly through L-Z 
and M*-Z diagnostics, improved models are crucial in pro- 
viding further insights into the evolution of local and distant 
galaxies. 

8. CONCLUSIONS 

With Spitzer IRAC imaging, we have derived total near- 
infrared luminosities for 27 nearby dwarf irregular galaxies at 
4.5 /im. The near-infrared L-Z relation is constructed using 
25 of these galaxies with secure distance and oxygen abun- 
dance measurements. While the slope of the subsequent [4.5] 
L-Z relation is slightly smaller than the corresponding opti- 



cal relations, the dispersion in the near-infrared L-Z relation 
is reduced compared to the optical L-Z relation. This result 
agrees with the expectation that variations in stellar mass- 
to-light ratios decrease from optical to near-infrared wave- 
lengths. We constructed a M»-Z relation, where the stellar 
masses were derived using color-based stellar mass-to-light 
ratios and corrected for the Chabrier IMF. The M»-Z relation 
for dwarf galaxies extends the SDSS M*~Z relation by an- 
other 2.5 decades. We find that the dispersion in the M*-Z 
relation is comparable over 5.5 decades in stellar mass. The 
dispersions in the L-Z and M*-Z relations as well as the large 
variation in effective yields are difficult to understand if galac- 
tic superwinds or outflows are responsible for low metallici- 
ties at low mass or luminosity. 
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TABLE 1 

Properties for 27 nearby (D < 5 Mpc) dwarf irregular galaxies. 



Nebular 







AORKEY 


%5] 


m-M 


Distance 




12+log(0/H) 


O/H 


B-[4.5] 


log M» 


Galaxy 


Source 


(G128) 


(mjy) 


(mag) 


Refs. 


(mag) 


(dex) 


Refs. 


(mag) 


(M ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


DDO 154 


S 




2.42 ±0.24 


27.53 ±0.10 


1 


-15.28±0.10 


7.67 ±0.05 


1,2 


1.88 


6.68 


DDO 165 


S 




9.57 ±0.96 


28.30±0.19 


2 


-17.62±0.19 






2.46 


7.73 


DDO 53 


S 




2.97 ±0.30 


27.76 ±0.15 


2 


-15.73±0.15 


7.62 ±0.20 


3 


2.36 


6.98 


GR8 a 


G128 


5054464 


2.95 ±0.30 


26.75 ±0.35 


3 


-14.72 ±0.35 


7.64 ± 0.04 


4 


2.52 


6.62 


Hoi 


S 




4.40 ± 0.44 


27.92 ±0.26 


1 


-16.32±0.26 


7.70 ±0.20 


5 


1.83 


7.08 


Ho II 


S 




42.9 ±4.3 


27.65±0.13 


1 


-18.53±0.13 


7.93 ±0.05 


6.7 


1.82 


7.96 


Ho IX 


s 




2.72 ±0.27 


27. 80 ± 0.08 b 




-15.75±0.10 


b 




2.07 


6.89 


IC 1613 


G128 


5051648 


91.53±9.2 


24.22 ±0.10 


4,5 


-15.92±0.10 


7.62 ±0.05 


8 


1.42 


6.82 


IC 2574 


S 




70.1 ±7.0 


28.02 ±0.22 


1 


-19.43 ±0.22 


8.15±0.11 


5.6 


2.09 


8.39 


IC5152 


G128 


5055232 


149 ±15 


26.58 ±0.1 8 


6 


-18.28zb0.18 


7.92 ±0.07 


8 


2.78 


8.11 


Leo A 


G128 


5052416 


8.65 ±0.87 


24.2 ±0.2 


7 


-13.34±0.20 


7.35 ±0.06 


4 


1.84 


5.89 


M81 dwB a 


S 




2.72 ±0.27 


28.62 ±0.10 


1 


-16.50±0.10 


7.98 ±0.22 


5 


1.96 


7.19 


NGC 1569 


G69 




220 ± 22 


26.71 ±0.60 


8 


-19.36 ±0.60 


8. 19 ±0.02 


9-11 


0.93 


8.07 


NGC 1705 


s 




17.4±1.7 


28.54 ±0.26 


9 


-18.44 ±0.26 


8.21 ±0.05 


12 


2.62 


8.13 


NGC 2366 


G69 




35.2±3.5 


27.52 ±0.28 


1,10 


-18.18±0.28 


7.91 ±0.08 


13,14 


2.18 


7.91 


NGC 3109 


G128 


5052928 


131 ± 13 


25.52±0.18 


11-13 


-17.61 ±0.10 


8.06 ±0.20 


7,15,16 


1.11 


7.41 


NGC 3738 


G69 




32.8 ±3.3 


28.45 ±0.25 


14 


-19.03 ±0.25 


8.23 ±0.01 


1 1 


2.45 


8.32 


NGC 4214 


G69 




175 ±18 


27.13 ±0.20 


15 


-19.53 ±0.20 


8.22 ±0.05 


17 


2.55 


8.55 


NGC 4449 


G69 




303 ± 30 


28. 12 ±0.27 


14 


-21.12±0.27 


8.31 ±0.07 


11,18 


2.98 


9.29 


NGC 5408 


s 




27.2 ±2.7 


28.41 ±0.17 


16 


-18.79±0.17 


7.98 ±0.01 


19,20 


2.29 


8.19 


NGC 55 


G128 


5056000 


785 ± 79 


25.85 ±0.20 


17-19 


-19.88 ±0.20 


8.05 ±0.10 


21-23 


1.58 


8.44 


NGC 6822 


S 




1440 ±140 


23.49 ±0.08 


20-22 


-18.14±0.10 


8.11 ±0.10 


24,25 


2.98 


8.10 


Peg DIG a 


G128 


5055744 


30.5 ±3.1 


24.90 ±0.10 


23 


-15.40±0.10 


7.93 ±0.13 


26 


2.90 


6.98 


Sextans A 


G128 


5053696 


8.39 ±0.84 


25.75 ±0.15 


24-26 


-14.85±0.15 


7.54 ±0.06 


1,27,28 


0.85 


6.24 


Sextans B 


G128 


5052672 


23.0±2.3 


25.63 ±0.04 


6,13,24,27 


-15.83±0.10 


7.53 ±0.05 


15,27-29 


1.73 


6.86 


UGC 6456 


G59 




2.25 ±0.23 


28.19 ±0.04 


13,28,29 


-15.87 ±0.10 


7.69 ±0.07 


11,14,30,31 


1.83 


6.90 


WLM 


G128 


5051136 


45.8 ±4.6 


24.83 ±0.08 


31-33 


-15.78±0.10 


7.83 ±0.06 


32 


1.88 


6.88 



NOTE. — Col. (1): Galaxy name (alphabetical order). Col. (2): Data source: G59 — GTO 59 (P.I. G. Rieke), G69 — GTO 69 (P.I. G. Fazio), G128 — 
GTO 128 (P.I. R. D. Gehrz), and S — SINGS (159; Kennicutt et al. 2003b). Col. (3): AOR Key pertaining to galaxies observed in program GTO 128. Col. (4): 
Flux density at 4.5 fim. Our measured flux densities for SINGS galaxies are consistent with those found in Dale et al. 12005); see also Fig.*3] Cols. (5) and 
(6): Distance modulus and references, respectively. Col. (7): Absolute magnitude at 4.5 fim. Cols. (8) and (9): Nebular (H II region) oxygen abundance and 
referen ces, respectively. Col. (10): Derived B— [4.5] colors in the present work, where B magnitudes are correct ed for Galactic e xtinctions; see Karachents ev et all 
1 2004). Col. (11): Derived stellar mass using model 8 in Table|3"| and adjusted downwards by 0.25 dex for the lChabrieJ <2oTj3T) IMF. 

RE FE RENCES. — Distance referen ce s: LIMakarova et al|19 98: 2. Karachentsev et al. 2002a; 3. Dohm-Palmer et al. 1998; 4. Cole e t alll999t 5.IDolphin et alj 
l200lt 6. IKarachentsey et alJI2002d 7. ITolstoy et aljll998t 8. lAngeretti et alj)2005t 9. ITosi et aUI2001t 10 . IThuan & Izotodl2005t lljMusella et alj|1997t 

12. IMinniti et alj|1999t~1 3. Mendez et al. 2002; 14. Karach entsev et alj|2003at 1 5. IDrozdovskv et al. 2002; 16. Karachentsev et al. 2002b; 17. Pritc het et all 
1 19871: 18. IKarachentsey et alJI2003bnTlTikhonov et alj|2005t 20. IGallart et alj[T99*a 21. IClementini et alj|2003t 22 IPietrzvnslri et alJI2004t 23. lApariciol 
WmL 24.TPiottoetalJll994t 2 5. ISakai et alj|1996t 26. IDolphin et alj|2003t 27. ISakai et alj|1997t 28. ILvnds et alj[T99*l 29. ISchulte-Ladbeck et alj|1998t 

3Q. IMinniti & Zii lstral 199*71 31. IDobhiri20001: 32. IReikuba et aH20001 

S pectroscopic oxygen a bun dance references: 1. Iv an Zee et al| 1 1997]: 2. [Kennicutt & Skillman 2001; 3. Skillman et al. 1989; 4. van Zee et al] 1200^; 
5. IMiller & Hodgdll996l: 6. IMaseeosa et aljlT99lt 7 . IH. Lee et aUI2003bt 8. IH. Lee et alj|2003at TlDevost et alj|1997t lO.TKobulmckv & Slcillmarin997t 
11. IMartirl 119971 12~Ih. Lee & Skillmarl 120041: 13. TG onzalez-Delgado et al. 1994; 14. Izotov et al. 1997t 15. IH. Lee et alj|2006bt 16. Richer & McCaU 
119951 : 17. Kobulnickv & Skillmanlll99'a 18. IBoker et alJI200U 19. IStasinska et aljll9861: 20. IMaseeosa et alJI 19941 21. ITalenflllTml 22. Iwebster & Smith! 
179831 23.frullmannetalJI2003t 24. IPeimbert et alj|2005t 25. IH. Lee et aUI2006at 26. ISkillman et alj|1997t 27. IKniazev et alJI2005r 28. lMagrini et 311120051 
29. IMoles et alll990t 3Q. lLvnds et al!1998l: 31. IHunter & Hoffmanll999t 32 Tk Lee et all2fJ05[ 



"Alternative names: GR 8 — DDO 155, UGC 80 91; M81 dwB — UGC 5423; Pega sus dwarf irregular (Peg DIG) — DDO 216, UGC 12613. 

b Ho IX: We have adopted the distance to M 81 i Ferrarese et al. 2000). Miller (1995) measured 12+log(0/H) ~ 8.0 ± 0.2 from a supernova remnant. 
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TABLE 2 

Comparison among selected luminosity-metallicity relations for nearby dwarf irregular galaxies and other distant 
star-forming galaxies. we show here a number of l-z relations from the literature; this is not intended to be a complete 
listing. we have not performed a homogeneous treatment of the data listed in this table. 





No. of 


Luminosity 


Wavelength 


L-Z fits 




Dispersion in 


Data source 


Galaxies 


Range 


(pm) 


a 


/? 




log(0/H) (dex) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


Local galaxies (D < 10 Mpc) 


This paper 


25 


-11 > M B >-18 


4.5 


5.78±0.21 


-0.122±0.012 


7.73 


0.12 


This paper 


25 


-11 >M B >-18 


0.4 (B) 


5.94 ±0.27 


-0.128 ±0.017 


7.99 


0.16 


Oliveira et al. 2006 


29 


-12 >M B >-18 


2.2 (K,) 


5.55 ±0.26 


-0.14 ±0.02 


7.79 


0.15 


van Zee et al. 2006 


33 


-11 >M B >-18 


0.4 (B) 


5.65±0.17 


-0.149 ±0.011 


8.03 


0.15 


H. Lee et al. 2003b 


22 


-11 >M B >-18 


0.4 (B) 


5.59 ±0.54 


-0.153 ±0.025 


8.04 


0.18 


Richer & McCall 1995 


21 


-11 >M B >-18 


0.4 (B) 


5.67 ±0.48 


-0.147 ±0.029 


8.02 


0.08 a 


Skillman et al. 1989 


20 


-11 >M B >-19 


0.4 (B) 


5.50 


-0.153 


7.95 


0.16 


Distant galaxies 


Rosenberg et al. 2006 


17 


-16 >M B >-18 


3.6 


4.85 


-0.16 


7.41 


0.16 


Salzer et al. 2005 


370 


-12 >M B >-22 


2.2 (K) 


3.92 ±0.09 


-0.212 ±0.003 


7.31 


0.24 


Savaalio et al. 2005 


79 


-18 >M S >-22 


0.4 (B) 


2.98 ±0.94 


-0.280 ±0.045 


7.46 




Kobulnickv & Kewlev 2004 


177 


-17 >M B >-22 


0.4 (B) 


4.90 


-0.193 


7.99 




J. Lee et al. 2004 


54 


-12 >M B >-20 


0.4 (B) 


5.37 ±0.46 


-0.159 ±0.029 


7.91 


0.26 


Tremonti et al. 2004 


> 53000 


-16 >M B >-22 


0.4 (B) b 


5.238±0.018 


-0.185 ±0.001 


8.20 


0.16 



NOTE. — Col. (1): Reference. Col. (2): Number of galaxies. Col. (3): Luminosity range. A number of the studies quote B,AB magnitudes, which are 
0.163 mag brighter than Johnson B magnitudes iFrei & Gunn 1994). Col. (4): Wavelength. Cols. (5) and (6): Fits to the respective relations are expressed as 
12+log(0/H) = a + bM\, where a and b are the intercept and the slope, respectively, and M\ is the absolute magnitude at wavelength A. Col. (7): «i6 represents 
the oxygen abundance determined from the specified relation at a fiducial galaxy luminosity equal to M\ = -16. Col. (8): Dispersion in oxygen abundance. 

11 Richer & McCall 1 1995) stated the fit parameters and resulting dispersion were effective only for galaxies brighter than M B = — 15, as the dispersion increased 
towards lower luminosities. However. IH. L ee et al. 1 2003b) showed that improved data (e.g., distances) reduced the dispersion at lower luminosities. 
b Tremonti et al. 1 2004) obtained galaxy luminosities in Sloan filter g (similar to B), and accounted for the appropriate transformations from g to B. 
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TABLE 3 

Stellar masses using models of stellar mass-to-light ratios from Tables 3 and 4 in IBell & de .TongH2001I) and the method 
described in §0 from the resulting derived stellar masses, we performed the resulting m»-z fit for each model. we list each 
m*-z fit as 12+log(0/h) = a + b loom*, and we list the dispersion in oxygen abundance for each fit. 







B- 


-V 


V-K 


B- 


-K 


M» 


-Zfits 


Dispersion in 


Model # 


Model 




b K 


UK 


b K 


a' 


V 


a 


b 


log(0/H) (dex) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Scaled Salpeter IMF, Z 


I 


Closed box 


-0.554 


0.540 


-0.740 


0.207 


-0.688 


0.150 


5.52 ±0.24 


0.307 ±0.031 


0.118 


2 


forfait 


-0.692 


0.699 


-0.926 


0.258 


-0.863 


0.188 


5.57 ±0.23 


0.304 ±0.030 


0.117 


3 


Outflow 


-0.534 


0.500 


-0.622 


0.157 


-0.601 


0.119 


5.49 ±0.24 


0.309 ±0.031 


0.118 


4 


Dynamical time 


-0.531 


0.476 


-0.687 


0.177 


-0.645 


0.129 


5.50 ±0.23 


0.309 ±0.030 


0.118 


5 


Formation epoch 


-0.694 


0.676 


-0.966 


0.270 


-0.888 


0.193 


5.58 ±0.23 


0.303 ±0.030 


0.117 


6 


Form, epoch: bursts 


-0.692 


0.652 


-1.087 


0.314 


-0.959 


0.212 


5.60 ±0.23 


0.302 ±0.030 


0.117 


7 


Cole et at. (2000) 


-0.654 


0.696 


-1.019 


0.301 


-0.909 


0.210 


5.58 ±0.23 


0.302 ±0.030 


0.117 








Different SPS 


a Models (IMFs); 0.4 Z Q 








8 


BC Salpeter 


-0.43 


0.60 


-1.16 


0.44 


-0.851 


0.254 


5.58 ±0.23 


0.298 ±0.030 


0.117 


9 


BC scaled Salpeter 


-0.59 


0.60 


-1.31 


0.44 


-1.005 


0.254 


5.62 ±0.23 


0.298 ±0.030 


0.117 


10 


BC modif. Salpeter 


-0.71 


0.60 


-1.45 


0.45 


-1.133 


0.257 


5.66 ±0.22 


0.297 ±0.029 


0.117 


11 


BC96 Scalo 


-0.65 


0.56 


-1.36 


0.43 


-1.052 


0.243 


5.64 ±0.23 


0.299 ±0.030 


0.117 


12 


KA Salpeter 


-0.38 


0.53 


-1.03 


0.39 


-0.754 


0.225 


5.55 ±0.24 


0.300 ±0.030 


0.117 


13 


Schulz et at, Salpeter 


-0.65 


0.77 


-2.50 


0.81 


-1.552 


0.395 


5.81 ±0.22 


0.284 ±0.030 


0.120 


14 


Pegase Salpeter 


-0.38 


0.59 


-1.15 


0.45 


-0.817 


0.255 


5.57 ±0.24 


0.298 ±0.030 


0.117 


15 


Pegase x = —1.85 


-0.07 


0.48 


-0.71 


0.37 


-0.431 


0.209 


5.45 ±0.24 


0.301 ±0.030 


0.117 


16 


Pegase x = —0.85 


-0.59 


0.71 


-1.25 


0.45 


-0.994 


0.275 


5.62 ±0.23 


0.296 ±0.030 


0.117 



NOTE. — Col. (1): Model number. Col. (2): Model type. BC - Bruzual & Chariot, KA - Kodama & Arimoto; see Bell & de Jond <200"il) for details. Cols. (3) 
and (4): Intercept and slope, respectively, for log (M*/Lk) as a linear function of B—V. Cols. (5) and (6): Intercept and slope, respectively, for log (M*ILk) as 
a linear function of V—K. Cols. (7) and (8): Derived intercept and slope, respectively, for log (M^ILk) as a linear function oiB-K. Cols. (9) and (10): Derived 
intercept and slope, respectively, for the subsequent M*-Z fits for the given model. Cols. (11): Dispersion in log(0/H). 

11 SPS models - stellar population synthesis models. 



TABLE 4 

Comparison among selected stellar mass-metallicity relations for nearby dwarf irregular galaxies and star-forming 
galaxies at distant epochs . fits to the respective relations are expressed as 12+log(0/h) = f(x), or functions of x = log m* . we 
have not homogenized measures of the stellar mass or metallicity for the other studies. this is not a complete listing, as we 

have selected recently published m*-z relations from the literature. 





No. of 




Luminosity 




M*-Z fits 


Dispersion in 


Data source 


Galaxies 


Redshift 


Range 




(Z = logM.) 


log(0/H) (dex) 


This paper 


25 


< 0.001 


-13>A%.5] > 


-21 


(5.54 ± 0.24) + (0.298 ± 0.030)X 


0.12 


Savaslio et al. 2005 


56 


~0.7 


-18>M fl ,,4 B > 


-22 


(4.06 ± 0.58) + (0.478 ± 0.058)Z 


0.20 


Tremonti et al. 2004 


> 53000 


~0.1 


-16 >M B >- 


22 


-1.492+ 1.847X-0.080X 2 


0.10 


Perez-Gonzalez et al. 2003 


163 


~ 0.026 


-19 >M K >- 


-26 


(6.02 ± 0.28) + (0.257 ± 0.028)Z 


0.21 
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TABLE 5 

Gas-to-stellar mass fractions and effective yields. 





r 2l 






Vrot 


fr H 1 H Y V 


( Tv Irm <;~ 1 A 
\J y Mil a j 


lop (M IM "1 


Ino A) cr 

lu & .Veff 


(km s~' ) 










C51 


DDO 154 


145 


2.00 


-1.26 


54 


DDO 165 


28.2 


0.55 




23 


DDO 53 


13.8 


0.770 


-2.51 


17 


GR8 


8.6 


0.464 


-2.75 


21 


Ho I 


38.9 


1.18 


-2.03 


18 


Ho II 


363 


1.17 


-1.82 


56 


Ho IX 


158 


1.93 




50 


IC 1613 


698 


1.26 


-2.04 


19 


IC 2574 


447 


0.97 


-1.90 


55 


IC 5152 


98.0 


0.02 


-2.84 


50 


Leo A 


68.4 


1.13 


-2.49 


4 


M81 dwB 


3.8 


0.34 


-2.52 


25 


NGC 1569 


74.1 


-0.01 


-2.59 


37 


NGC 1705 


16.6 


0.02 


-2.55 


61 


NGC 2366 


295 


1.07 


-1.93 


45 


NGC 3109 


1 1 10 


1.39 


-1.47 


51 


NGC 3738 


21.9 


-0.09 


-2.61 


50 


NGC 4214 


324 


0.32 


-2.30 


43 


NGC 4449 


794 


0.36 


-2.18 


84 


NGC 5408 


64.6 


0.49 


-2.37 


33 


NGC 55 


2679 


1.15 


-1.72 


78 


NGC 6822 


2399 


0.17 


-2.54 


55 


Peg DIG 


22.8 


-0.36 


-3.07 


8 


Sextans A 


208.8 


1.83 


-1.56 


49 


Sextans B 


102.4 


0.92 


-2.45 


22 


UGC 6456 


14.1 


1.03 


-2.19 


7 


WLM 


299.8 


1.03 


-2.05 


26 



NOTE. — Col. (1): Galaxy name. Col. (2): 21-cm flux densities i Karachentsev et al. 20 03). Col. (3): Gas-to-stellar mass ratio. Col. (4): Logarithm of the 
effective yield. Col. (5): H I rotational velocity corrected for inclination and turbulent motions I Karachentsev et al. 2004). 
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FIG. 1. — Postage-stamp images at 4.5 fj,m of dwarf galaxies from GTO 128. North and East are at the top and left, respectively, in each frame. All galaxies 
are displayed with the same isophotal scaling to emphasize the large range in surface brightness of the dwarf galaxies. The fields of view shown are: GR 8 — 
5' X 5', IC 1613 — 10' X 10', IC 5152 — 5' X 5', Leo A — 5' X 5', NGC 55 — 15' X 10', NGC 3109 — 10' X 5', Peg DIG — 5' X 5', Sextans A — 
5' X 5', Sextans B — 5' X 5', and WLM — 5' X 15'. 



FIG. 2. — Example of galaxy photometry for Sextans B. North and East are to the top and left, respectively, in the 4.5 fim image. The central elliptical aperture 
surrounds the galaxy, and the seven circular apertures sample the background flux. The galaxy aperture is aligned roughly to the same orientation as the optical 
emission seen on images from the Digitized Sky Survey. 
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FIG. 3. — Comparison of total 4.5 ^m fluxes for selected dwarf galaxies from the SINGS sample: fluxes from Dale et al. 1 2005) versus measured fluxes from 
the present work. The dotted line marks the line of equal fluxes. 



14 



H. Lee et al. 




-12 -14 -16 -18 -20 -22 

M [4.5] ( ma §) 

FIG. 4. — (a) B absolute magnitude versus [4.5] absolute magnitude. Lines of constant B-[4.5] color between and 4 are also plotted, (b) B-[4.5] color versus 
[4.5] absolute magnitude. The Pearson correlation coefficient for this color-magnitude relation is r = -0.24. 




FIG. 5. — Nebular oxygen abundances versus (a) absolute magnitude in B, and (b) absolute magnitude at 4.5 fim. The solid lines represent fits to the data, 
whose fit parameters are given in each panel. The Pearson coefficients are r = -0.80 in B, and r = -0.89 in [4.5]. 
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FIG. 6. — (a) Residuals in oxygen abundance (from B fit) versus B absolute magnitude, (b) Residuals in oxygen abundance (from [4.5] fit) versus [4.5] absolute 
magnitude. Both panels are plotted to the same vertical scale, and a short-dashed horizontal line in each panel indicates zero residual. The dotted horizontal lines 
in each panel indicate the 1<t scatter. The scatter in the NIR L-Z relation is reduced compared to that found in the B L-Z relation. 
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FIG. 7. — Panel (a): Residuals in oxygen abundance from the B L-Z fit versus B-[4.5] color. Panel (b): Residuals in oxygen abundance from the [4.5] L-Z 
fit versus i?— [4.5] color. Panels (a) and (b) are plotted to the same vertical scale. Panel (c): Residuals in B absolute magnitude from the B L-Z fit versus £-[4.5] 
color. Panel (d): Residuals in [4.5] absolute magnitude from the [4.5] L-Z fit versus B— [4.5] color. Panels (c) and (d) are plotted to the same vertical scale. 
Residuals at optical wavelengths correlate (r = +0.67) more strongly with B-[4.5] color than residuals at NIR wavelengths (r = +0.28). 
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12+log(0/H) = 
(5.65 1 0.23) + 
(0.398 ± 0.030) log M, 
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FIG. 8. — Panel (a): oxygen abundance versus stellar ma ss. Stellar masses have been derived using the color-based stellar mass-to-light ratios from model 
8 in Tablel3l and have been corrected to the Chabriei 1 2003) IMF. The fit to the data is shown as the solid line; the Pearson coefficient is r = +0.90. Panel (b): 
residuals in oxygen abundance from our present M t -Z relation against stellar mass. The Pearson coefficient is r = —0.08. 




6 7 8 9 10 11 12 

log M„ (M ) 

FIG. 9. — Oxygen abundance versus stellar mass. The plot has been expanded to include more massive galaxies. Data and the least-squares fit to the present 
data from Fig. |8lar e shown; the fit is shown as a thin black solid line and is arbitrarily extrapolated to 1O IO M0. The polynomial fit to over 53000 galaxies 
from the SDSS i Tremonti et al. 2004) is shown as a heavy red solid line; the la curves are also shown. Their fit is effective over the mass range between 8.5 < 
log (M*/Mq) < 11.5. In the "overlap region" (8.5 < log (M./Mq) < 9.5), SDSS galaxies are more metal-rich by PS 0.2-0.3 dex compared to dwarf galaxies in 
the present sample. An extrapolation of the SDSS M„-Z relation down to log (M*/Mq) = 7.5 is shown as a blue short-dashed line. 
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FIG. 10. — Panels (a) and (b): oxygen abundances are plotted against gas-to-stellar mass ratio. In panel (a), the solid line represents the fit to the data; the 
Pearson coefficient is r = —0.58. In panel (b), we have plotted outflow variations of the closed-box model, where the outflow rate is a constant multiple of the star- 
formation rate. The purple and blue curves describe models with yo = 0.01 i Tremonti et al. 2004) and yo = 2 X 10~ 3 (e.g., H. Lee et al. 2003b; van Zee & Havnes 
2006), respectively. For each set of curves, the outflow rates from right to left are 0, 5, 10, 25, and 50 times the star-formation rate, respectively. In panel (c), 
B— [4.5] colors are plotted against gas-to-stellar mass ratio. The solid line represents the fit, and the Pearson coefficient is r = —0.83. The fit excludes DDO 165, 
Ho IX, and NGC 1569; the latter is the outlier with a very blue color for its relatively low gas-to-stellar mass ratio. 
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Fig. 11. — Panel (a): effective yield against total baryonic mass in gas and stars; cf. also lGarn gt3(20 i O i 2 [). Filled circ les represent dwarf galaxies 
in the present work . Crosses represent th e median of the SDSS data in mass bins of 0. 1 dex jTrernorltietal] 2004). The empirical fit is shown 
as a thin solid line jTremonti et a l. 2004, their E quation 6). The asymptotic yield (y e ff = 0.0104) is represented by the s hort-dashed line. The 
solar oxygen mass fraction l Asplund et al. 2004; Melendez 2004) and the effective yield for the nearby spiral galaxy M5 1 I Bresoli n'et all2004T) 
are represented by dotted and long-dashed lines, respectivel y. These value s are comparable to the range of effective yields for dwarf galaxies 
(e.g., log y c ff — -2.7 to -2.1. IH. Lee et all2003b! ; lvan Z ee & Havnes 2006). Panel (b): effective yield versus gas-to-stellar mass ratio. The fit 
sh own is repr esented by a thin solid line. The trend here is similar to the one observed for a sample of nearby UGC dwarf galaxies described 
bv lvan Zee & Havnes! ( 2006 ). The horizontal lines are the same as in the panel above. 
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FIG. 12. — Panel (a): stellar mass versus H I rotation velocity. The solid line represents the fit to the data with Pearson coefficient r — +0.67. Panel (b): gas- 
to-stellar mass ratio versus H I rotatio n velocity. Panel (c) : effec tive yield versus H I rotation velocity. H I rotation velocities in TablelTlhave been corrected for 
inclination and turbulent motions; see Kar ache ntsev et al. 1 2004). The Pearson correlation coefficient between rotation velocity and effective yield is r = +0.40. 
The horizontal lines are the same as in Fig. II II 
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